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We test the response to a threading magnetic field of rings of 5-site C4„-symmetric repulsive Hubbard 
clusters connected by weak intercell links; each 5-site unit has the topology of a Cu04 cluster and 
a repulsive interaction is included on every site. In a numerical study of the three-unit ring with 
8 particles, we take advantage of a novel exact-diagonalization technique which can be generally 
applied to many-fermion problems. For O-O hopping we find Superconducting Flux Quantization 
(SFQ), but for purely Cu-Cu links bound pair propagation is hindered by symmetry. The results 
agree with W = Q pairing theory. 
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The repulsive Hubbard ring in the presence of a magnetic 
flux has been studied by several authors [1] [2] [3] [4]. An 
anomalous Aharonov-Bohm (AB) effect [5] with ground-state 
energy oscillations versus flux (j) having a period shorter than 
the foundamental one = hc/e was reported, but no sign 
of superconductivity was found. The rings that we consider 
in this paper have nodes which consist of 5-site units with a 
Cu04 topology. In the following we shall refer to the central 
site of each node as Cu and to the four external sites as O 
just to distinguish their position in the unit cell, even if we 
acknowlegde that the present model may be too idealized to 
have much direct relevance to high-Tc cuprates. The 5-site 
cluster will be also called CUO4. Below we show numerical 
solutions of such a model that clearly show superconducting 
pair hopping if the total number of particles is 2|A| +2p where 
jAj is the number of units and < p < |Aj; in particular, once 
a magnetic field is switched on into the ring, SFQ is unam- 
biguously observed. We note that in our model no super- 
conducting response is obtained with less than 2 particles per 
5-site unit; the Zhang-Rice picture [6] for the two-dimensional 
d — p model does not represent superconducting pairs and for 
the present small system we need to explore a scenario with 
a slightly larger density of particles. By replacing CUO4 by 
larger units one can model other ranges of filling fraction. 

The repulsive Hubbard Hamiltonian has two-body singlet 
eigenstates without double occupation [7] [8] [9] [10] called 
W = Q pairs. Such solutions are also allowed in the fully 
symmetric clusters C. In the many-body ground state these 
pairs get dressed and bound, and this is signaled by Ac(A^) < 
where Ac(A^) = Ef\N) + e'^°\N - 2) - 2E^°\N - 1); 
E^^ (N) is the interacting ground state energy of the cluster 
C with N particles. By means of a non-perturbative canonical 
transformation [11] [12], it can also be shown that Ac (A'') < 
is due to an attractive effective interaction and at weak 
coupling I Ac (A'') I is just the binding energy of the pair. The 
extension of the theory to the full plane was also put forth in 
Ref. [11]. The dv symmetric 5-site cluster is the smallest one 
where the W = pairing mechanism works. We have already 
described W = pairing in great detail as a function of the 
one-body and interaction parameters on all sites; the study 
was extended to larger clusters too [9] [13]. 

We are using CUO4 as the unit just for the sake of simplic- 



ity, but the W — mechanism produces bound pairs at differ- 
ent fillings for larger clusters [10] too. In order to simplify the 
analysis, here we neglect the 0-0 hopping term within each 
unit and the only nonvanishing hopping matrix elements are 
those between an O site and the central Cu site; they are all 
equal to t. In the total Hamiltonian 

Htot ~ Ho + Hr, (1) 

Ho describes the units and H-r the hopping between first- 
neighbor units. 

|A 

Ho = ^[i^(dLPa,i<T +piAad'a^) + 
a — l ia 

(^(^i't'^S+E'^ST^a)] (2) 

i 

where p], ^^ is the creation operator onto the O i — 1, .., 4 of 
the a-th unit and so on. The point symmetry group of Ho 
includes S^^^ , with |A| the number of nodes; in this report we 
consider |A| = 2 and |A| = 3. 

We take Hr invariant under the Si subgroup of S^^^ , al- 
though a square symmetry would be enough to preserve the 
Acu04 (4) < property. We analyzed alternative models for 
this term. 

First,we studied hopping between corresponding O sites, 
i.e. connecting the i-th O site of the a-th unit to move towards 
the i-th O site of the /3-th unit with hopping integral Tc/s = 

-ffr = E E [Tc,fjpl_i^Pl3,i,T + h.c] . (3) 
a — /3il ia 

In numerical work we assumed small (tc/jI ^ jAcu04(4)|. 
This is necessary because the clusters that we can diagonalize 
explicitly are small and a large |r| with a flux of the order of 
a fluxon would produce a strong perturbation on the ground 
state multiplet. The same flux in a large system would indeed 
be perturbative even for |r| ~ t. 

For N = 2|A| and Tap = 0, the unique ground state consists 
of 2 particles in each CUO4 unit. We assume a total number of 
N = 2|A|-|-2p particles. When U/t is such that Acu04 (4) < 0, 
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each of the p added pairs prefers to He on a single Cu04 and 
for A'' = 2|A| + 2p the unperturbed ground state is 2''x (''^') 
times degenerate (since ^£ has dimension 2). 

We exactly diagonalize the |A| = 2 and |A| = 3 ring Hamil- 
tonian. The numerical work is made easy by the powerful 
Spin-Disentangled technique, which we briefly introduced re- 
cently [14], but deserves a fuller illustration. To the best of 
our knowledge, it was not invented earlier, which is somewhat 
surprising, being a very general method. 

We let Mf + Mj^ = A'' where is the number of particles 
of spin u; {l^aa)} is a real orthonormal basis, that is, each 
vector is a homogeneous polynomial in the and d) of degree 
Mct acting on the vacuum. We write the ground state wave 
function in the form 



|*> =^L„;3|<A„T>®l<A/3i> 



(4) 



a/3 



which shows how the | and | configurations are entangled. 
The particles of one spin are treated as the "bath" for those of 
the opposite spin: this form also enters the proof of a famous 
theorem by Lieb [151. In Eq.(4) Laj3 is a m| x my rectangular 
matrix with mCT=(^'). We let denote the kinetic energy 

rua X rua square matrix of -fftot in the basis {\4>a<j)}, and Ns'^'^ 
the spin-a occupation number matrix at site s in the same 
basis {Ns"^ is a symmetric matrix since the |0c«t)'s are real). 
Then, L is acted upon by the Hamiltonian iftot according to 
the rule 



(5) 



In particular for M| = Mi (Sz = sector) it holds = Ki 
and A^i^-* = Ns^\ Thus, the action of H is obtained in a spin- 
disentangled way. The generality of the method is not spoiled 
by the fact that it is fastest in the Sz ~ sector, because it 
is useful provided that the spins are not totally lined up; on 
the other hand, Sz = can always be assumed, as long as the 
Hamiltonian is SU(2) invariant. 

For illustration, consider the Hubbard model with two sites 
a and b and two electrons (H2 molecule) each in the (j>a or <j)b 
orbital. The intersite hopping is t and the on-site repulsion 
U. In the standard method, one sets up basis vectors for the 
Sz = sector 



1^1) = l(*aT) «> \<Pal), 



\^2) = \M<g>\<pbl), 
IV4) = \M 8> \<t>bl), 



One then looks for eigenstates (three singlets and one triplet) 



(6) 



of the Hamiltonian 
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(7) 



Insted of working with 4x4 matrices, we can cope with 2x2 
by the spin-disentangled method using the form in Eq. (4) with 



L = 



N, 



1 




Ka = 



N, 



t 

t 
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Using Eq.(5), one finds 

H»2W =Y.Y1 (HhALMM ® IM) (8) 



cx = a,b j3 = a,b 



with 



Ulpl + tilp2 + 1p3) 



Utpi + t{tp2 + Vs) 



(9) 



The reader can readily verify that this is the same as applying 
-H"h2 in the form of Eq.(7) to the standard wave function in 
Eq.(6) and then casting the result in the form of Eq.(4). Since 
we can apply Hyi2 wc can also diagonalize it. The advantage 
of working with 2x2 rather than 4x4 matrices for the H2 toy 
model is ridiculous, but it grows with the size of the problem 
and in the |A| = 3 ring case it is spectacular. In the 5^=0 
sector for |A|=3 the size of the problem is 1863225 and the 
storage of the Hamiltonian matrix requires much space; by 
this device, we can work with matrices whose dimensions is 
the square root of those of the Hilbert space: 1365 x 1365 
matrices solve the 1863225 x 1863225 problem, and are not 
even required to be sparse. 

Here we have implemented this method for the Hubbard 
Hamiltonian. We emphasize, however, that this approach will 
be generally useful for the many-fermion problem, even with 
a realistic Coulomb interaction, which can be suitably dis- 
cretized. Starting from a trial wave function of the form in 
Eq.(4) we can avoid computing the Hamiltonian matrix, since 
its operation is given by Eq.(5). Each new application of the 
Hamiltonian takes us to a new Lanczos site and we can pro- 
ceed by generating a Lanczos chain. To this end we need 
to orthogonalize to the previous sites by the scalar product 
given by (^'i|^'2) = Tr(LjL2). In this way we put the Hamil- 
tonian matrix in a tri-diagonal form. This method is well 
suited since we are mainly interested in the low-lying part of 
the spectrum. A severe numerical instability sets in when the 
chain exceedes a few tens of sites, i. e. well before the Lanczos 
method converges. Therefore we use repeated two-site chains 
alternated with moderate-size ones. 

In the basis of the sites (of the original cluster) the oc- 
cupation matrices A^i"^^ are diagonal with elements equal to 
or 1, simplifying the calculation of the interaction term. 
Moreover, in choosing the trial wave function for small r we 
take full advantage from our knowledge of the 5*4 irrep of the 
T = ground state. This speeds the calculation by a factor 
of the order of 2 or 3 compared to a random starting state 
(or even more, if U is large). Typically, starting from a t = 
ground state for the three-unit ring, 24 short Lanczos chains 
were enough to obtain a roughly correct energy and a 20- site 
chain achieved an accurate eigenvalue and an already stabi- 
lized eigenvector. In limiting cases when the results could 
be checked against analytic ones, using double precision rou- 
tines an accuracy better than 12 significant digits was readily 
obtained. 
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The two-unit ring (|A| = 2) does not allow to insert a 
flux; the reason is that each unit is at the left and at the 
right of the other, and one cannot tell which is the clock- 
wise motion. Results on the spectrum of such a system will 
be reported in a more comprehensive paper. Here, for the 
three-unit ring (jA| — 3) we focus on the case y < 34 (which 
grants Acu04(4) < 0) and p = 1 (total number of particles 
2|A|-F2p = 8). 

The inter-unit hopping r between the O sites breaks the 
symmetry group Csu Sf into Civ S4 for real r; in a mag- 
netic field (complex t), this further breaks into C31S1S4. More 
explicitly, let us label the ground state multiplet components 
by the crystal momentum 2nhk/3 where k is an integer. Real 
r separates = from the k = 1 and k = 2 subspace of 
C3 , but cannot split A; = 1 and 2 because they belong to 
the degenerate irrep of C31,; this degeneracy is resolved by 
complex r , when we insert a magnetic flux </!> by r = |T|e'*, 

The persistent diamagnetic currents carried by bound pairs 
screen the magnetic flux. We calculated the expectation value 
of the total current operator [16] 

i,a.,a 

(10) 

as a function of the flux and versus k; this current operator 
yields a gauge invariant average /. We note incidentally that 
by expanding / in Eq.(lO) in powers of (j) near (f> = one may 
identify paramagnetic and diamagnetic contributions with the 
zeroth and the first order terms respectively [17], although 
such a splitting is not gauge invariant. 

The numerical results for / are reported in Fig.(l); they 
also convey direct evidence of the SFQ since according to 
the Hellmann-Feynman theorem the current is proportional 
to the flux derivative of the ground-state energy. At = 
the ground state has fc = 0. Near (j> = the system gener- 
ates a diamagnetic current which screens the threaded mag- 
netic field; = is a local minimum of the ground state en- 
ergy, which grows quadratically in (j> (diamagnetic behaviour) . 
However, when <j> exceedes a critical value ~ <^o/4, a symme- 
try change of the ground state occurs due to a level crossing 
to = 2. This corresponds to a sharp discontinuity of the 
current which suddenly changes sign; thus, for stronger flux 
the current enhances the external field until, ai (j> ~ <^o/2 
the current vanishes again, signaling the energy minimum at 
half fluxon. This anti-screening of the fleld beyond the level 
crossing is interesting but the paradox is easily understood. 
Indeed, like at = 0, the eigenfuctions at </!> = (j)o/2 may be 
choosen real, since Ht{4> ~ 4>o/2) is obtained from Ht{4' ~ 0) 
by reversing the sign of four O-O r bonds connecting two 
nearest neighbours units. Thus, near (f>Q/2 the magnetic flux 
can still be considered as a small perturbation, but the un- 
perturbed real inter-unit Hamiltonian is Hr{<i> = 4>o/2); the 
current normally screens this perturbation. The vanishing of 
the current at <^ = (j>o/2, when the ground state energy is 
in a new minimum belonging to the k — 2 subspace, also 
marks the restoring of the time reversal invariance, like in 
the BCS theory [18]. Indeed, at (p — (j)o/2 the symmetry 
group is Csv <8> <S'4 where Csv is isomorphous to Csv (re- 
flections a are replaced by ag, where <? is a suitable gauge 



transformation). This feature was also found in other ge- 
ometries [11] [13]. Finally, a level crossing at a critical value 
~ 3(^0 /4 takes to the k = 1 ground state and to the mini- 
mum at = (/>o; at this point the flux can be gauged away 
taking us back to i^i = 0, fc = 0. We numerically verified 
that Ek=2{4>) = -Efc=2(0o - 0), Ek=o{(t)) = Ek=i{<Po - (j>) and 
Ek=i{<t>) ~ Ek=o{4>o — (t>), where Eki(f>) is the ground state 
energy in the k sector. Thus, the dressed W = pair screens 
the vector potential like one particle with an effective charge 
e* = 2e does. At both minima of E^'^\(j}) we have computed 
A3_unit(8) ~ —10~^t. Here, the half-integer AB effect is ac- 
tually SFQ. From Fig.(l) we see that the maximum value of 
the diamagnetic current is of the order of 1 10 nano Am- 
pere ii t = 1 eV and the ratio ^^^^ « e\r\/h near = 0. 
The present results are in line with the expectations of the 
analytic theory that we have presented elsewhere [14]. As 
foreseen, this pattern disappears when U/t —> 0. 
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FIG. 1. Total current for the three-Cu04 ring, as a function of 
the magnetic flux. Here U = 5 t, \t\ = 0.0014. The current is in 
units of e|r|/h. The thick line marks the ground state current. 

Next, we report the effects of an inter-unit hopping rcu 
between Cu sites only; this does not break the S]^^ symmetry 
and therefore its consequences on pair propagation and SFQ 
are drastically different. In order to study the propagation of 
a bound pair we again assume the total number of particles 
2|A| +2 — 8. The full system threaded by the fiux has a 
C3 ® Sl symmetry because the O sites are not involved in 
the inter-unit Hamiltonian. We computed the ground state 
energy versus flux for the case Acu04(4) < 0; rcu produces 
much smaller effects than r for I rcu I ^ |Acu04(4)|. Therefore 
we are free to take the inter-unit hopping as large as |rcul = 
O.lt without causing a drastic change; the dependence of the 
ground state energy is still weak, see Fig. (2). At = 0, the 
correction [19] due to |tcu| to ground state energy is only 
~ 10~'^t. However, the most spectacular change is that the 
system now behaves like a paramagnet. 
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the m X m problem by means of a suitable disentanglement 
of the up and down spin configurations. 
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FIG. 2. Ground state energy E of the three-unit ring in 
units of t, as a function of the concatenated magnetic flux. E 
is fc-independent (see text). Here U = 5 t, \tcu\ = O.lt. 



The reason of this unusual behavior is the local symmetry, 
which forbids any fiux-induced splitting of the three k levels; 
the effective mass is infinite and the W — pair, however 
bound, is localized. Indeed the S4 label of each CUO4 unit 
is a good quantum number. No SFQ is observed because the 
screening of the magnetic field by the bound pair is forbidden. 
The small correction to the ground state energy comes from 
a second-order process. Suppose we prepare the system e.g. 
in an unperturbed state |4, 2, 2), that is , with 4 particles in 
the first Cu04 and 2 in the others; the bound pair is localized 
on the first cluster. The evolution of the wave packet involves 
virtual states |3, 3, 2) and |3, 2, 3) which are obtained when a 
particle jumps to the nearby clusters. This particle must be 
in totalsymmetric state, because the irreps of each unit must 
be conserved. However, such a process occurs with a tiny 
amplitude because the energy misfit is severe. This is partic- 
ularly clear at weak coupling, when we can speak in terms of 
orbitals, and the lowest-energy Ai particle of the first cluster 
must hop to antibonding Ai orbitals of the nearby clusters; 
the amplitude of this process is further reduced by the over- 
lap of these orbitals with the localized Cu one. Moreover, 
such virtual processes are insensitive to the flux. Any (p de- 
pendence arises from third-order corrections in which the Ai 
particle virtually goes around the trip clockwise or anticlock- 
wise. In the ground state, of course, it chooses the wise in 
such a way to gain energy from the magnetic field. This is 
why a paramagnetic dependence on the flux is seen in Fig. (2) 
and the correction goes like — c/i^ at small (j). This is interesting 
because it shows how the local symmetry can hinder the tun- 
neling of bound pairs carrying conserved quantum numbers; 
SFQ is not a necessary consequence of superconductivity if 
the pairs are not totalsymmetric. 

In conclusion, we have found that Hubbard-like repulsive 
graphs hosting W = 2-body states can be used to model 
bound pairs and their propagation. Thus, SFQ may be found 
in purely repulsive Id Hubbard models only if the nodes are 
represented by a non-trivial basis. For rings of Cu04-units 
and weak O-O links we find a half-integer AB effect which is 
unambiguosly interpreted as SFQ. Next, we reported a coun- 
terexample, namely, the case of Cu-Cu links, when pairing is 
not leading to SFQ due to the large effective mass of the parti- 
cles. By manipulating matrices whose dimension is the square 
root of the overall size of the Hilbert space we were able to ex- 
actly diagonalize the Hamiltonian: y'm x y'm matrices solve 
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